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Abstract. The causal perturbation expansion defines the Dirac sea in the presence 
of a time-dependent external field. It yields an operator whose image generalizes the 
vacuum solutions of negative energy and thus gives a canonical splitting of the so- 
lution space into two subspaces. After giving a self-contained introduction to the 
ideas and techniques, we show that this operator is in general not idempotent. We 
modify the standard construction by a rescaling procedure giving a projector on the 
generalized negative-energy subspace. The resulting rescaled causal perturbation 
expansion uniquely defines the fermionic projector in terms of a series of distribu- 
tional solutions of the Dirac equation. The technical core of the paper is to work 
out the combinatorics of the expansion in detail. It is also shown that the fermionic 
projector with interaction can be obtained from the free projector by a unitary trans- 
formation. We finally analyze the consequences of the rescaling procedure on the 
light-cone expansion. 
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1. Introduction 

Shortly after the formulation of the Dirac equation j2j, it was noticed that this 
equation has solutions of negative energy, which have no obvious physical interpretation 
and lead to conceptual and mathematical difficulties. Dirac suggested to solve this 
problem by assuming that in the physical vacuum all states of negative energy are 
occupied by electrons forming the so-called Dirac sea O Hj . Since this many-particle 
state is homogeneous and isotropic, it should not be accessible to measurements. Due 
to the Pauli exclusion principle, additional particles must occupy states of positive 
energy, thus being observable as electrons. Moreover, the concept of the Dirac sea led 
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to the prediction of anti-particles. Namely, by taking out particles of negative energy, 
one can generate "holes" in the Dirac sea, which are observable as positrons. 

Today, Dirac's intuitive concept of a "sea of interacting particles" is often not taken 
literally. In perturbative quantum field theory, the problem of the negative-energy 
solutions is bypassed by a formal replacement and re-interpretation of the creation 
and annihilation operators of the negative-energy states of the free Dirac field, giving 
rise to a positive definite Dirac Hamiltonian on the fermionic Fock space. In the 
subsequent perturbation expansion in terms of Feynman diagrams, the Dirac sea no 
longer appears. This procedure allows to compute the S'-matrix in a scattering process 
and gives rise to the loop corrections, in excellent agreement with the high-precision 
tests of quantum electrodynamics. 

One shortcoming of the perturbative approach is that the particle interpretation 
of a quantum state gets lost for intermediate times. This problem becomes apparent 
already in the presence of a time- dependent external field. Namely, as first observed 
by Fierz and Scharf [5], the Fock representation must be adapted to the external field 
as measured by a local observer. Thus the Fock representation becomes time and 
observer dependent, implying that also the distinction between particles and anti- 
particles loses its invariant meaning. The basic problem can be understood already in 
the one-particle picture: Mathematically, the distinction between particles and anti- 
particles corresponds to a splitting of the solution space of the Dirac equation into two 
subspaces. In the vacuum, or more generally in the presence of a static external field 
"B(x), in the Dirac equation 



The separation constant u, having the interpretation as the energy of the state, gives 
a natural splitting of the solution space into solutions of positive and negative energy. 
The Dirac sea can be introduced by occupying all states of negative energy. However, 
if the external field is time dependent, 



the separation ansatz no longer works, corresponding to the fact that the energy of 
the Dirac states is no longer conserved. Hence the concept of positive and negative 
energy solutions breaks down, and the natural splitting of the solution space seems to 
get lost. 

Another shortcoming of the standard reinterpretation of the free Dirac states of 
negative energy is that this procedure leads to inconsistencies when the interaction is 
taken into account on a non-perturbative level. For example, in [131 [T4] the vacuum 
state is constructed for a system of Dirac particles with electrostatic interaction in the 
Bogoliubov-Dirac-Fock approximation. In simple terms, the analysis shows that the 
interaction "mixes" the states in such a way that it becomes impossible to distinguish 
between the particle states and the states of the Dirac sea. Thus the only way to obtain 
a well-defined mathematical setting is to take into account all the states forming the 
Dirac sea, with a suitable ultraviolet regularization. 

The framework of the fermionic projector is an approach to formulate quantum 
field theory in such a way that the above-mentioned problems disappear. Out of 
all the states of the Dirac sea we build up the so-called fermionic projector, which 



(iVfy + ®(x) - m) *(ar) = 
one can separate the time dependence with the plane-wave ansatz 

V(t,x) = e~ iwt ij(x) . 
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puts Dirac's idea of a "sea of interacting particles" on a rigorous mathematical basis. 
The fermionic projector gives a global (i.e. observer-independent) interpretation of 
particles and anti-particles even in the time-dependent setting at intermediate times. 
The interaction is described by an action principle, which can be formulated without 
referring to the causal or topological structure of the underlying space-time, thus giving 
a possible approach for physics on the Planck scale (see [12] for a review and more 
references). When analyzed in Minkowski space, this action principle yields all the 
Feynman diagrams of perturbative quantum field theory, but also gives rise to other 
surprising higher order corrections [10]. The foundations of this approach have been 
worked out in the book [9]. The connection to the fermionic Fock space formalism and 
to second quantized bosonic fields is elaborated in [TT] . 

The construction of the fermionic projector is based on the observation first made 
in [6] that for the time-dependent Dirac equation (jl.ip there still is a natural splitting 
of the solution space into two subspaces, if one works instead of the sign of the energy 
with the underlying causal structure. To explain the basic idea (for details see Section [2] 
below), we first note that the image of the operator with integral kernel 

J*T (*, y)=j^L + m ) S tf _ m 2) e -Hx-y) (L2) 

coincides precisely with all the negative-energy solutions of the free Dirac equation. 
This operator can be decomposed as 

where the operator k m is causal (in the sense that its kernel is supported inside the 
light cone), and the operator p m can be understood as the absolute value of the op- 
erator k m . In the case with general interaction, one can extend k m uniquely to an 
operator k m using the causality property. Making sense of the absolute value, one can 
also generalize p m to an operator p m . Introducing in analogy to (11.30 the operator 

Pr = l(p m -~km), (1-4) 

the image of this operator describes the Dirac sea in the case with interaction. 

The aim of the present paper is to clarify the normalization of the states of the 
generalized Dirac sea (|1.4|) . More precisely, the operator (jl.2p describing the free Dirac 
sea has the property that it is idempotent if a ^-normalization in the mass parameter 
is used, 

This idempotence property plays an important role in the framework of the fermionic 
projector. However, as we shall see, the interacting Dirac sea (jl.4p as defined in [6] is 
in general not idempotent in this sense. Our goal is to modify the normalization of the 
states using a rescaling procedure such as to arrange (jl.5p . This issue of normalizing 
the states of the interacting Dirac sea can be regarded as a problem of functional 
analysis. To see the analogy, if the Dirac operator were a self-adjoint operator on a 
Hilbert space, we could interpret the product dm as the operator- valued spectral 
measure of the Dirac operator, composed by a projector on the generalized negative- 
energy solutions of the Dirac equation. Unfortunately, the Dirac operator is only 
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symmetric with respect to the indefinite inner product 



<^|<I>> 



/ 



V(x)$(x) d A 



x 



(1.6) 



making it impossible to use spectral theory in Hilbert spaces. This is the reason why 
we must rely on perturbative techniques and work with formal power series expan- 
sions. The main technical task is to work out the combinatorics of the perturbation 
expansions in detail. The interesting point is that the details of these expansions have 
a correspondence to general results known from functional analysis in Hilbert spaces. 
In particular, we recover the polar decomposition, the resolvent identity, Stone's for- 
mula and the functional calculus from our perturbation expansions (see (|2.7p . (|2.19p . 
Remark 14.21 and (|2.25p ). We also relate the free and interacting operators by an oper- 
ator U which is unitary with respect to the indefinite inner product (jl.6p (see (|5,ip ). 

The main result of the paper is the derivation of a unique perturbation expansion 
for the fermionic projector which satisfies (jl.5p (see Theorem 14. ip . The summands of 
this expansion can be regarded as Feynman tree diagrams, as they also appear in the 
standard perturbation expansion in the presence of an external field. However, our 
perturbation expansion is different in that the usual freedom in choosing the Green's 
functions (like working with the advanced or retarded Green's functions or the Feyn- 
man propagator) is removed. The expansion becomes unique by combining causality 
with suitable normalization conditions for the states of the fermionic projector. 

The paper is organized as follows. In Section^ we give a self-contained introduction 
to the construction of the generalized Dirac sea (II. 4p in terms of a formal perturbation 
expansion in S. In Section [3l we explain the rescaling procedure for the states of 
the interacting Dirac sea, thus obtaining a unique idempotent operator in terms of 
a formal perturbation expansion, the so-called rescaled causal perturbation expansion. 
The main technical task is to elaborate the combinatorics of the different expansions in 
detail; this will be carried out in Section HI As explained in Remark 14. 2\ the rescaling 
formally reproduces results from spectral theory and functional analysis in the setting 
of perturbation expansions. An interesting consequence is the existence of unitary 
transformations between the free and the generalized Dirac seas. We prove this fact 
in Section [5] by deriving and analyzing equations for the perturbation flow. Finally, 
in Section [6] it is shown that the rescaling procedure has no influence on the residual 
argument and the light-cone expansion as worked out in [8]. But it does change the 
form of the so-called high-energy contribution. 



In this section we give a self-contained review of the causal perturbation expansion 
as developed in [6] (see also Chapter 2]). The rescaling procedure will then be 
explained in Section [3j We always assume that the mass is a positive parameter, 



2. The Causal Perturbation Expansion 



m > . 



We decompose the Fourier integral (II. 2p as follows 
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where 



Pm{x, y) = 1 4\± (0 + m ) 5 ^ 2 - ™ 2 ) e ~ tq{x ~ y) 

J ( J (2-2) 



(2k)* 

Here d = q^^, e is the sign function, and the function in (jl.2p is the Heaviside 
function. For the signature of the Minkowski inner product we use the convention 

(H ). All these Fourier integrals are well-defined tempered distributions. The 

splitting (12. ID gives rise to the decomposition (II. 3h of the corresponding operators. 

The decomposition (|2.ip reveals the following connection to causality. The Dirac 
equation is causal in the sense that information propagates at most with the speed of 
light. This is reflected in a support property of the advanced and retarded Green's 
functions, which we denote by and s^, respectively. They have the Fourier repre- 
sentation 

sl(x,y)= f^L * + m 



my ,yJ J (2tt) 4 q 2 -m 2 - ieq " 

s a (x v) = f $ + m e -M*-v) 

a m \x,V) J (2?r)4 q2 _ m 2 + i£qo 

where e > is a regularization parameter, and it is understood implicitly that one 
should take the limit e \ in the distributional sense. Computing the integrals with 
residues, one readily verifies that supp(s^(x, .)) C J% , where 

= { y with (y - x) 2 > and y° > x } 

denotes the future light cone centered at x (similarly, supp(s^(x, .)) C J£). Taking 
the difference of the two expressions and using the identity 

= i (zrKz ~ zr 1 -) . (2-3) 



(2.4) 



2m' \ x — is x + is 

one finds that 



km — ^ . ( S ™ S m) 



This relation shows that k m is a causal operator in the sense that supp(fe m (x, .)) C J x , 
where J x = J% U j£ is the light cone centered at x. 

We point out that the operator p m is not causal in the above sense. To see this, we 
decompose p m in analogy to (|2.4|) as 

Pm = ~2^l^ Sm ~ s m) ' (2-5) 
where the operators s ± have the Fourier representation 

J (z7r)* q z — m z =p is 

The operator s~ is known in the literature as the Feynman propagator, character- 
ized by the condition that positive-energy solutions propagate forward in time while 
negative-energy solutions propagate backwards in time. An explicit calculation of the 
Fourier integral (|2.6|) in terms of Bessel functions shows that s+ and s~ as well as 
their difference do not vanish outside the light cone and are thus not causal. 
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For the subsequent constructions it is important to observe that the operator p m 
can be obtained from k m as follows. Defining the absolute value of a diagonalizable 
matrix A as the unique positive semi-definite matrix \A\ with A 2 = \A\ 2 , we find that 
\e(k°)(fc + m)\ = (ft + m). Since the operators p m and k m are diagonal in momentum 
space, taking the absolute value pointwise in momentum space can be understood 
formally as taking the absolute value of the corresponding operator acting on the 
Dirac wave functions, 



Pm — \km\ • (^■'O 



In the remainder of this section we shall generalize the relations (|2.4p and (|2.7p to 
the case with general interaction (jl.lj) ; in Section [3] we will then develop a method for 
generalizing the definition of the fermionic projector (jl.3j) . 

Using the causal support property, the advanced and retarded Green's functions s m 
and s m are uniquely defined even in the case with interaction (jl.ip . This could be 
done non-perturbatively using the theory of symmetric hyperbolic systems |15j . For 
our purpose, it is sufficient to give the unique perturbation series 



S~m ~ / / ( s m-^) s m ' s m — / v ( s m 23) s m . (2.1 
n=0 n=0 

Here the operator products involving the potential S are defined as follows, 

(s m ^s m ){x,y) := d*z s m (x, z)"B(z)s m (z,y), (21 



It is straightforward to verify that the perturbation series (]2.8p indeed satisfy the 
defining relations for the advanced Green's function 

(i$ x + 'B(x) - m) s^(x, y) = 5 4 (x - y) , supp s^(x, .) C J x 

(and similarly for the retarded Green's function; for details see [9j §2.2]). 
Moreover, we define the operators 

1 v A i 

s m = 2^ Sm ^ Srn ^ = 2^ Sm Sm ^ , (2.10) 

F m (Q,n) = \ Pm > UeQ for Q C N, n € N (2.11) 

I k m : n £ V 

and the series of operator products 

oo oo oo 

b m = ^(s^r , b m = £(-S Sm )"B , b> = £(-S* m ) n . (2.12) 

n=0 n=0 n=0 

Lemma 2.1. The following identities hold: 

PmPm' = k m k m > = 5(m - m') p m (2.13) 

Pm k m i = k m p m , = 5(m- m') k m (2.14) 
PP 

Pm Sm' — Sm' Pm — 7 Pm (2.15) 

m — ml 

PP 

km Sm' — Sm' k m — 7 k m (2.16) 

m — m 

PP 

s m s m ' = 7 (s m - s m >) + ir 2 5(m - m!) p m , (2.17) 
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where the principle value distribution is given by — := h[(x + ie)^ 1 + (x — ie)^ 1 ]. 
Proof. Calculating pointwise in momentum space, we obtain 

Pm(q)Pm'(q) = % 2 - m 2 )($ + m)5(q 2 - m 2 ){f + m) 

= 5(m 2 — m' 2 )5(q 2 — m 2 ) (q 2 + (m + m')$ + mm'^j 

= S(m — m')5(q 2 — m 2 )(m 2 + (m + m')d + mm') 

2m r 

= 5(m — m')5(q 2 — m 2 )2m(m + d) = 5{m — m')p m (q) . 

2m p 

This gives the first part of (|2.13p . The second part of this formula as well as formula 
(|2.14p are obtained analogously. The formulas (|2.15p and (|2.16p are computed in the 
following manner: 

2 o ( 4 + m ' d + m' 

2p m (q) s m i{q) = 5(q - m )(d + m) -= -z : h 



q 2 — m' 2 — ieqo q 2 — m' 2 + ieqo 
5(q 2 - m 2 )(q 2 + (m + m')d + mm') ( \ : h 



q 2 — m' 2 — ieqo Q 2 — m' 2 + i^Qo 



6(q 2 — m 2 )(m 2 + (m + m')d + mm') ( — ^ -~ : h 



5{q 2 - m 2 ){d + m) ( 
PP 

■jPm(q) ■ 



m 2 — m' 2 — ieqo m2 ~ m ' 2 + i£Qo 
(m + m') j (m + m') 



(m + m')(m — m') — ieqo ( m + m')(m — m') + ieqo 



m — m 



The derivation of (|2,17p is a bit more difficult. By (|2.4|) and (|2.1U|) . we have 



Sm — Sjn Mrfem — ^k m • 



Thus we can express the product s m (q) s m i{q) in two ways, namely as 
s m (q) s m ,{q) =( s m(<?) - mk m (q))(s y m , (q) -mk m ,(q)) 

= s m(?) s m'(?) - <7r ( k m'(q)— - : + k m (q) : J 

\ m — m — ieqo m — m — ieqo J 

- ir 2 5(m - m')p m (q) , 

or alternatively as 

s m (q) s m >(q) =( s m(q) + i ^k m (q))(s^ n ,{q) + mk m >{q)) 

=»m(«) s^ n ,(q)+iTT (k m ,(q)— 5— + k m (q) 

\ m — m + ieqo m — m + ieqo 

- ir 2 5(m - m')p m (q) . 
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Adding these two formulas yields 

2s m (q)s m >{q) - (s^{q) s y m , (q) + s^{q) s^, (q)) + 2vr 2 5(m - m')p m (q) 

= iitk m i (q) 



1 1 



+iirk m (q) 



m + zeqo ml — m — isqo 
1 1 



m — m + ieqo m — m — ieqo 



== mk m i (q)e(-qo) 2iri 5(m' — m) + i7rk m (q)e(—qo) 2iri 5(m — m) 
= -27r 2 5(m' - m)(-p m >(q)) - 2i?8{m - m!){-p m {q)) , 



where in the last step we used the definitions of p m and k m . We thus obtain 

1 
2 

It remains to derive the relations 



*„> *,„' - - (sm sjjj, + s^,) + ir 2 5(m - m!)p m . ( 2 -18) 



V V V V A A A A 

s m s m' = ~ ~j( S m ~ s m>) an( ^ s m s m' = ~~ Zj( S m ~ s m') > (2-19) 

m — ml m — ml 

which can be regarded as "resolvent identities" for the causal Green's functions. It 
suffices to consider the case of the advanced Green's function. Clearly, the operators 
on the right side of (12. 19ft satisfy the support condition supp((s)^ — s^ n ,)(x, .)) C J% , 
and from 

s m s m'( x ^y)=J d A zs y m (x,z) s y m ,{z,y) 

we see that the operators on the left side of (|2.19p satisfy this support condition as 
well. Moreover, the calculations 



and 



(i$ x - m) s^ n s w m , (x, y) = s y m , (x, y) 



PP 

(i$ x - m) 7 (s^ - s m ,)(x,y) 

m — ml 

PP 

7 [6(x -y)- {m! - m)s^ n ,(x,y) - 5(x - y)] = s^,(x,y) 



m — m 



show that both sides of (|2. 19|) satisfy the same inhomogeneous Dirac equation. Hence 
their difference is a distributional solution of the homogeneous Dirac equation vanishing 
outside jv. The uniqueness of the solution of the Cauchy problem for hyperbolic PDEs 
yields that this difference vanishes identically. This shows ()2. 19[) and thus finishes the 
proof of (I2T7D. □ 

We remark that the last summand in (|2.17p was by mistake omitted in [9]. For 
this reason, we shall now rederive and correct the formulas of the causal perturbation 
expansion in detail. 

Corollary 2.2. Let C S {p m , k m } and C £ {p m ', k m /} and , as in (|2.12p . Then 
the following calculation rule holds: 

Cb>b< f C = CC' + 5{m - m') tt 2 C b mPm b m C. (2.20) 
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Proof. Using the calculation rules of the previous lemma, we obtain 

/ 1 \ PP 

C[ y"('Bs m ) l (s m/ 'B) n - l )C > = CsmrtC' + CBsmC' = (CBC - CBC) =0. 

V — ' / m — m 

1=0 

The same method also applies to higher order. We then get a telescopic sum, but the 
last summand in (|2.17f) gives additional contributions. More precisely, for any n > 2, 



n 

c(Y,('% S m) l (Smi'B) n - l )c' = 
1=0 

n-l 

= C{Bs m ) n C + C0w3) n C" + C [ Y,(Zsm) l (s m >'Z) n - 1 } ( ' ' 

i=i 

PP 



C — L -^—f [-{Ks^-H + S^B)"- 1 ] C 
m — m' 

n-l 



/ pp \ 

+ C ^{'BsrJ-H 7 (s m - s m + 5(m - m')-n 2 p m ^{s^T' 1 ' 1 C 

7 c [-(s^)"- 1 ^ + nsm^r- 1 ] c 



1=1 

PP 



m — rrv 

pp n-l n-2 



+ 1 c(Y d ^3 m ) l ^s m ,y^ l -H - Y^eBs m ) l (<Bs ml ) n - l - l K)C' 

1=1 1=0 
n-l 

+ 5(m - m')vr 2 C ^(Ss m ) i - 1 Bp fn B(s m /B) n -'- 1 C" 
l=i 

n-l 

-- 5(m - mV C ^(B 5m ) i - 1 Bp m B(5 m ,B) ri - i - 1 C" 
l=i 

n-2 

-- 5(m - mV C Y,C%s m ) l 'Bp m 'B{s m ,'B) n ~ l - 2 C l . 



1=0 

We thus obtain 



cb>b<, c = cj2(-^ m r J2(-sm>v) n 'c 

n=0 n'=0 

oo n 

n=0 1=0 

oo /n—2 \ 

= CC' + S(m - m')7r 2 ^(-l) n C7 ^(S Sm ) / S Pm S( Sm ,S) n - z - 2 C 



n=2 \J=0 
oo / n 



CC' + 5{m - mV £(-l)»C ^(S Sm ) z S Pm S( Sm ,S)^ ] C 

n=0 \l=0 

CC' + 5{m — m')ir 2 C b m p m b m C' '. 
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We next generalize the operator k m to the case of non-static external potentials. 
The natural way to do this is to define 

' = 2~7r7^ m _ ' (?-^) 

giving a rigorous non-perturbative definition of k m . The formal perturbation expansion 
of k m is given in the following proposition. 

Proposition 2.3. The relation (|2.2ip uniquely determines the perturbation expansion 
for k m . We have the formula 

CO 

~k m = ^2(-i7r) 2f3 b<k m (b m k m )^b> . (2.22) 

/3=0 

Proof. An explicit calculation shows that 

{%$ + S - m) 6< = . (2.23) 

As all operator products in (12. 22ft have a factor at the left, the series in (|2.22j) is a 
solution of the Dirac equation. 
From fl23J) and (|2TTU|) . we have 

— ^ m ^fom ) — s m i7rk m . (2.24) 

We substitute the series f|2.8j) into (|2.2ip . insert (|2.24p and expand. A reordering of 
the resulting sum gives the claim. The details of the reordering process can be found 
in 0. □ 

In Appendix [A] we give a compilation of the perturbation series of k m and of other 
perturbation series which will appear in what follows, up to third order in H. 

Next, we want to generalize the operator p m to the case with interaction. Note that 
for the construction of k m , it was crucial that k m was a causal operator. However, Prn 
does not have this property, and it is not at all obvious how to generalize p m using the 
decomposition (|2.5p . Namely, as mentioned after (jl.ip . in the time-dependent setting 
the energy is not a conserved quantity, and thus the frequency conditions characterizing 
the Feynman propagator (see after ([2. 6ft ) have no meaningful generalization. Instead, 
we want to exploit equation (12. 7p to generalize the operator p m , i.e. 

Pm " y \f&L ■ (2.25) 

This will be done using a Taylor series for the square root function, defining p m by a 
perturbation series. 

Using Corollary 2.2 and Proposition 2.3, the operator k m k m ' is calculated to be 
k m k m ' = 5{m — m) 



oo 



^ v m)-^ • -™b<(kmb m ) 2l3l (p m + TT 2 k m b m p m b m k m \(b m k m ) W2 b^ l . (2.26) 



In the following calculations, all operator products can be computed using the rules 
of Lemma |2. II and Corollary 12.21 This always yields a factor 5{m — m 1 ), and thus for 
notational simplicity we can omit these factors 5{m — m') and consider all expressions 
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at the Scinie value of TTi. Furthermore, leaving out the subscripts i m ', our calculation 
rules can be written in the compact form 

p 2 = k 2 =p (2.27) 

pk = kp = k (2.28) 

b > b < = p + TT 2 pbpbp . (2.29) 

In order to keep the combinatorics simple, it is convenient to first treat the sum- 
mand ir 2 pbpbp in (|2.29p by rearranging the perturbation series of k. 

Definition 2.4. We define the series of operator products 

b< ■= b< (p + ^pbpbp)- 1 ' 2 = b<(p + V(-l) n (2n ~ l )" (TT 2 pbpbp) n ) (2.30) 

V 2 n n\ ) 

n>l 

b> := (p + ^pbpbp)- 1 ' 2 b> = ( p + V (-1)" (2n - (ir 2 pbpbp) n ) b> (2.31) 

V 2 n n\ / 

n>l 

p + A:={p + ^pbpbpfl 2 = p + ^(-1)" +1 {2n ~ 2 T ^ 2 pbpbp) n . (2.32) 

n>l 

Manipulating formal power series, like for example 

b>&< = (p + ■K 2 pbpbp)~ 1/2 b > b < (p + ir 2 pbpbp)~ 1/2 

= (p + iT 2 pbpbp)~ l l 2 (p + TT 2 pbpbp)(p + iP'pbpbp)^! 2 = p , 
we find the simpler calculation rules 

b > b < =p, (2.33) 
b < (p + A) = b < , (2.34) 
(p + A)b> = b> , (2.35) 
(p + A) 2 = p + i?pbpbp . (2.36) 

Thus we can rewrite k and & 2 as 

oo 

k = ^2(-in) 2l3 b < (p + A)k(bk) 2f3 (p + A)b > , (2.37) 

oo 

P = ^ {-in) 2Wl+f32) b< (p + A) k{bk) Wl (p + vr 2 ' pbpbp)k{bk) W2 (p + 4)&> . (2.38) 

/3i,/3 2 =0 

We are now ready to derive the perturbation expansion for p. 

Theorem 2.5. TTie relation (|2.25p uniquely determines the perturbation expansion for 
p. We have the formula 

P = b< U + X>* n j b> (2.39) 
wii/i i/ie combinatorial factors 

cn = (-ir +l (2ra ~ 3)!! (2-40) 
v ; n!2 n 
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and the operator 

oo 

X = -p+ ^2 {-^) 2Wl+l32) (p + A)k(bk) 2 ^(p + A) 2 k(bk) 2 ^{p + A). (2.41) 

Proof. Using the calculation rule (|2,33p . all the intermediate factors b < and b > will 
always drop out of our calculations. Therefore, for ease in notation we can simply omit 
all factors b < and b > . Following (|2.25p . we are thus looking for a positive operator p 
being a powers series in 23, such that 

CO 

p 2 = k 2 ™ Yl (-iK) m+M (p + A) k(bk) 2 ^ {p + n 2 pbpbp)k(bk) 2 ^ (p + A) . 

Using again the operator X defined in (|2.4ip . we have 

p 2 =p + X. (2.42) 

The operator p is idempotent and acts as the identity on X. Thus we can take the 
square root of p + X with a formal Taylor expansion, 

P = VP~TX = p + E(- 1 )" +1 ^?^ n > ( 2 - 43 ) 

n>l 

which uniquely defines p as a positive operator. Reinserting the factors b < , ir into 
d233|), we obtain the formula d239|). □ 

Note that, due to (|2.30p and (|2.23p . the series (|2.39p is a solution of the Dirac equa- 
tion (jl.ip . Furthermore, rewriting (|2.4ip as 

X = 2A + A 2 + ir 2 (kbpbk + Akbpbk + kbpbkA + AkbpbkA) 

+ Yj {-iTT) 2Wl+l32) {p + A)k(bk) 2f3l (p + Tr 2 pbpbp)k(bk) 2 ^{p + A), (2.44) 

(/3i,/W(0,0) 

one sees that X only involves operator products with at least two powers of b (and thus 
also at least two powers of 23). Consequently, X n is of order 0(23 2n ), so that (|2.39[) is 
indeed a well-defined power series in 23. 

The following Lemma shows that under suitable regularity and decay assumptions 
on the potential 23, all operator products appearing in (|2.22p and (|2.39p are well-defined 
and finite. 

Lemma 2.6. Let n £ N and Cq,Ci, . . . ,C n G {k,p,s}. If the external potential is 
smooth and decays so fast at infinity that the functions 23(x), XiB(x) and XiXjH(x) 
are integrable, then the operator product 

(C„23C n _ 1 23...23C )(x,y) (2.45) 
is a well-defined tempered distribution on M 4 x M 4 . 

For the proof, one derives Schwarz norm bounds by combining estimates in position 
and in momentum space. The details can be found in {9j Lemma 2.2.2]. 
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3. The Rescaling Procedure 



Introducing in analogy to (|2.ip the operator 



t:=i(p-fc), (3.1) 

the key feature of the operator t is that it specifies via its range a certain subspace of 
the space of solutions of the Dirac equation (jl.ip . Since in the vacuum case, the range 
of t is precisely the space of all negative-energy solutions, we refer to the range of t as 
the generalized negative- energy solutions. 

In 0[9], the interacting fermionic projector is defined by (|3.ip . However, the op- 
erator t is not a projection operator in the sense of (II, 5j) . This can already be seen 
by considering the perturbation expansions to second order. Namely, from (12.37p 
and (1^391) . 

t = -b < (p-k + kbkbk + -kbpbk + -pbpbp - -pbkbk - -kbkbp) 6> + 0(6 4 ) 
2 V 2 2 2 2 / 

Taking the square and using the calculation rules (j2.33H2.36]) . we obtain 

t 2 = —b < (^p — k + —kbkbk + —kbpbk + —pbpbp — pbkbk 

— kbkbp + -pbkbp — -kbpbp — -pbpbk^j Ir* + 0(6 4 ) . 

Obviously, these expansions for t and t 2 do not coincide. 

Our strategy for resolving this problem is to rescale the states in the range of t with 
the following general procedure. Since the fermionic projector P should have the same 
range as the operator t, we take the ansatz 

P = tYt 

with a "rescaling operator" Y which should be invertible and should commute with t. 
We define 

R = it and Z = kk=pp. (3.2) 

The simplest idea would be to choose U Y = i? -1 ", but as the range of t are only the 
generalized negative-energy solutions, the operator R will in general not be invertible. 
On the other hand, the operator k maps onto the whole space of solutions of the Dirac 
equation, and thus the formal definition "Y = Z" 1 " does the trick. We will now 
make this definition precise and show that it indeed yields a projection operator. As 
explained on page LTO] we again omit all factors b < and b > , knowing that they can be 
reinserted at the very end. According to (|2.42|) . Z is given by Z = p + X. As p is 
idempotent and acts as the identity on the operator X, we can define Y by a Neumann 
series, 

Y =p + Y J {-V l X l . (3.3) 

l>\ 

Then 

ZY=(p + X){p + Y J {-Xf)=p, (3.4) 

l>i 

showing that Z and X are inverse in our context where p acts as the identity on all 
other operators. 
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We are now ready to define the fermionic projector in the interacting case (jl.ip and 
to prove that we indeed obtain a projection operator. 

Definition 3.1. The fermionic projector P is defined by 

P:=tYt. (3.5) 
with t according to (|3.ip and Y given by (|3.3|) . 
Proposition 3.2. The fermionic projector is idempotent, 

P 2 = P , 

(where, as explained on paae\T(K we leave out the factor 5(m — m') and compute for a 
fixed mass parameter m). 

Proof. The formulas (|2.37|) and (|2,41[) show that the operators k and X commute, 
k{p + X) = (p + X)k 

oo 

= {-i7r) 2{ ^ 2+P3 \p + A)k{bk) 2 ^{p + A) 2 k(bk) 2 ^(p + A) 2 k{b^ . 

Pi A, 

Noting that the operators p and Y are formed as a sum of powers of X (cf. (|2,43|) 
and (|3.3|) ). it follows that these operators both commute with each other and with k, 

\p,Y] = \p,~k] = [Y,k] = 0. 

Hence 

_ _ 1 ~ ~ 1 ~~ ~~ 1 

tYt = -(p - k)Y(p - k) = -{ppY + kkY - pYk - kYp) = -(p - pYk) , (3.6) 

and thus 

PP = (tYt){tYt) = -(p - pYk){p - pYk) 

1 _ _ _ i 

= -(p - 2pYk + pYkpYk) = j(p- ^pYk + (ppY)(kkY)) 

= -(p - 2pYk + p)= tYt = P . n 



4. Derivation of the Perturbation Expansion for the Fermionic 

Projector 

The definition of the fermionic projector f)3.5[) has the disadvantage that it involves 
several hierarchies of infinite sums. In this section we shall simplify the perturbation 
expansion of the fermionic projector by bringing it into a form with the same structure 
as the perturbation expansion (|2.22p for k involving only two hierarchies: P should be 
an infinite sum of operator products of the form b < FibF2b...bF n b > with Fj G {p, k}, 
where the operators b < , b, b > are defined by the perturbation series (|2.12p . Again, the 
leading orders of the expansions derived in this section can be found in Appendix lAl 

Theorem 4.1. The fermionic projector P m defined in (|3.5[) can be written as 

Pm = \(V™-K S )i (4-1) 
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where the operators p™ s and k™ s are given by 

oo 

Pm = &m Y.^) W PmOmP m ) W b> (4.2) 

(3=0 
oo r 

c s = b< J2 E c ^ p)°^ p) 6 - • ( 4 - 3 ) 

r=0 p=0 

Here we have used the definitions 



and 



G m (r, P )= (-l)' 7(r * P * Q) F m (Q, 1) b m F m (Q, 2) b m b m F m {Q, 2r + 1) , (4.5) 

QeT(2r+X) 
#Q=2p 

where V(2r + 1) denotes the power set of {l,...,2r + 1}, the operators F(Q,n) are 
defined in (12.11j) and 

a(r,p,Q) = l+ x - C 46 ) 

xe{l,...,2r+l}\g 

Here the superscript "res" indicates that the operators p Tes and k res have been rescaled 
using the procedure of Section [3l Alternatively, the superscript "res" can be regarded 
as an abbreviation of the word "residual," in view of the "residual argument" used in 
the light-cone expansion of the fermionic projector (cf. Section EJ). 



Proof of Theorem\4--l\ We begin with formula (13.60 . 



P = - b< (p - kpY) b> . (4.7) 

Using that p = (p + X) 1 / 2 and Y = (p + X)^ 1 , we obtain 

pY = (p + X)-2 = p + e n X n = £ tnX n 

n>l n>0 

with the combinatorial factors 

(Or, _ 1 VI 

= ^ —FnT 1 ' n > 1 • (48) 
For the last identity, we have also set eo = 1 and X° = p. Inserting into equation 
(14. 7jl . we obtain 



P = I b< p b> - - b< J2 e n ~kX n b> (4.9) 

n>0 

It is convenient to give the first series in (|4.9p a name, 



p™ s := b< p b> = b< E(i7r) 2 ^(5p) 2/3 b> . (4.10) 
One can show in a way completely analogous to Proposition 12.31 that 
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where s + and s~ are the power series that arise formally analogous to (|2.8p . 

oo oo 

^m = ^""^(~ s m-^) s m ) = ^]( — s m^) s m ■ (4-12) 

n=0 n=0 

The remaining task is to simplify the formula for the second series in (|4.9p . which 
we now denote by £; res , i.e. 

£re S ._l< J^dnkX 71 ^ , (4.13) 
n>0 

This will be done in Lemma 14.31 concluding the proof of the theorem. □ 

Remark 4.2 (Stone's formula and polar decomposition). In this remark we ex- 
plain formal analogies of our perturbation expansion with formulas known from func- 
tional analysis. First, for a selfadjoint operator I? in a Hilbert space, one can compute 
the spectral projector T m using Stone's formula (see e.g. [161 Theorem VII. 13]), 

y m = -L {{V - (m + ie))- 1 — (T> — (m — ie))" 1 ) . (4.14) 

Formally applying this identity to the vacuum Dirac operator T> = i0, the resolvent is 
a multiplication operator in momentum space; namely, 

\ -l / \ — 1 d + m 

V-(m±ie)) (q) = (t-(m±ie)) = q2 _ {m±l£ y 

Comparing with (|2.5p and (|2.6p . we find that the spectral projector coincides with the 
operator p m with integral kernel (|2.2p , 

1 
2~vri 



Tm — Pm — n __. (s m S m ) 



In the case of the interacting Dirac equation V = i$ + 23, the resolvent can be 
computed formally with a Neumann series, 



(p-(m±ie)) 1 = J](-s ± S) n s d 



— a U j f 

n>0 

This coincides precisely with the perturbation series (I4.12p . and thus we can write 
Stone's formula (|4.14p as 

— fir - .s- ) - r ,Vfl 
2m 



— r,_-(^m *m) ~ Pm ■ 



with as given by (|4.1ip . This consideration shows that the operator p Tes defined in 
(I4.10P should be regarded as the spectral projector of the interacting Dirac equation. 
This is consistent with the equations 

p res k res = k res = k res p Tes and p res p res =p res = k Tes k res , (4.15) 

which can be computed explicitly from (I4.10P and (I4.13p . 

Next, a bounded linear operator K in a Hilbert space admits a polar decomposition 
of the form 

K = U\K\, 
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where \K\ = y/ K*K is the absolute value of K and U is an isometry from the or- 
thogonal complement of ker(K) to the closure of the image of K (see e.g. [1 6|. Theo- 
rem VI. 10]). From (|4.7p . we see that 

k res = kpY = kp- 1 = k\k\~ l . 

We therefore regard fc res as the isometry transforming \k\ into k and having the same 
range and kernel as k. Since p is a positive operator, i.e. p = \p\, and p res is the spectral 
projector and thus acts as the identity on the range of p, one can also regard p Tes as 
the isometry transforming \p\ into p and having the same range and kernel as p. 

It remains to provide the lemma quoted in the proof of Theorem 14.11 

Lemma 4.3. The operator k res defined by (I4.13P can be written as 

oo r 

k ies = b<^2^2c(r,p)G(r,p)b> (4.16) 

r=0 p=0 



with the coefficients 



2r r(r- P + |) 



c{r,p)=^ \ " '» (4.17) 
T(-p + i)H 



and the operator products 



G(r,p)= Yl (-l) a{r ' p ' Q) F(Q,l)bF(Q,2)b---bF(Q,2r + l) . (4.18) 



QeV(2r+l) 
#Q=2p 

Proof. We introduce the abbreviation 
U :=p + X 

cm 



{-iTT) 2Wl+l32) (p + A)k(bk) 2 ^(p + A) 2 k{bk) 2l32 {p + A) . 

/3i,/3 2 =0 

Since p acts as the identity on U, we can calculate the n-th power of X to be 

= (-p +ur=j2 (") (- P ) n - l u i = p (") (-I)-' u i , 



where the 1-th power of U is given by U° = p and 

t ,-i»)-- , k(bk)-->i»+ Ar 
if / > 0. Moreover, from (|2.37p we have 



U l = {p + A) (J^ {-mf P k{bkf p (p + A) 2 y * ^ Q (-ivr) 2Q A:(6A:) 2Q (p + A) 



k = Yi-i^ip + A)k(bk) 2/3 (p + A) . 

13=0 
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We thus obtain 



kX n = ^( U \-l) n - l kU l 
1=0 ^ ' 



ivr) 2/3 fc(6fc) 2/3 (p + ^)' 



2/ 



x V (-i7r) 2Q fc(6fe) 2a (p + ^) . 
z — /« 

Inserting into (|4.13p and using the calculation rules (|2.34p - (|2.36p . we obtain the fol- 
lowing formula for k res : 



n=0 Z=0 ^ 7 



2/ 



x 6< ^(-ivr) 2/3 A:(6A:) 2/3 (p + ir 2 pbpbp) ^{-m) 2a k{bk) 2a b> . (4.19) 



<f3=0 



a=0 



Remarking that in (|4.19p one gets one factor tt for each order of b, for notational 
convenience we shall omit the factors tt during the proof and reinsert them at the very 
end. Note moreover that for a given n, the /-sum in (|4.19p equals kX n and is therefore 
of order 0(b 2n ) (as X is of order 0(6 2 ), cf. (|2.44p ). Thus for every r G N, the truncated 
sum 



n=0 1=0 



1) 



n—l 
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xb< \J2(-lfk(bk) 2p (j) + pbpb P )j ^{-l) a k(bk) 2a b> (4.20) 

coincides with k res up to terms of order 0(6 2r+2 ). For this reason, we can calculate 
^.res f rom (|4 20p to every order in b. Interchanging the n- and Z-sums (which are both 
finite), we can carry out the n-sum, which involves only coefficients but no operator 
products, to obtain 



k 1 



£ 

1=0 



,n=l 



\n—l 



b< ...b> + 0(6 



2r+2^ 



2/ 



Y,^b< \ Y,{-ifk{bk) 2 ^{p+pbpbp)\ ^2(-irk(bk) 2a b> + o(b 2r+2 ) 



1=0 

with the coefficients 



fl,r = J2e n (^J* 



-I) 



(2r + l)!! 
l\ 2 r (2l + l)(r-l)\' 



The combinatorics of the above operator products is analyzed in Lemma 14.41 below. 
Using the results of this Lemma, we obtain that for all r 6 N, the contributions to A: res 
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involving 2r factors of b are given by 

r I 



/=0 p=0 v r ' QG-P(2r+1) 

#Q=2p 

x F(Q, 1)6F(Q, 2)6 •• • 6F(Q, 2r + 1) 6> 

^EE^f't^W,^ ( 4 - 21 ) 

(=0p=0 \ 1 P / 

(where we used the abbreviation (14.180 ). Summing over all orders in 6, we conclude 
that 

oo r I / i 7 \ 

^ S = ^EEE Ar ( " T_ " P J G(r, p) 6> . (4.22) 

r=0 Z=0 p=0 ^ ^ ' 

Since the operator products G do not involve the index /, we may interchange the sums 
over I and p and perform the Z-sum, 



r 



+ Z-/A r(r-p + |) 



l-W T(i-p)r! 



[Y,(-^fk(bk)^(p + pbpbp)] 2l ^2(-irk(bk) 2a , (4.23) 



where we used the formula [U eq. (6.1.12)] to write the double factorial in terms of the 
gamma function with half-integer argument. Inserting back the factors it, we obtain 
the result. □ 

We finally prove the following combinatorial lemma. 

Lemma 4.4. For any p,l,r £ No with p < r, the following statements hold. 

(i) In the expression 

\{-lfk{bkfP(p + pbpbp) % 



only terms with an even number of factors p appear. 

(ii) In the series (|4.23|) . each term of the order b 2r of the form 
C l bC 2 b...C2rbC 2 r + i , CiE{p,k} (4.24) 

appears exactly ( r ^~ p ) times, where 2p denotes the number of factors p in the 

product. Here we adopt the convention (",) =0 if n' < 0. 
Moreover, each term of the form (|4.24p appearing in (|4.23p has the sign 

(-1)^sj x where J = {i € {1, 2r + 1} with Q = k} . 
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Proof, (i) We proceed by induction in I. For I = 0, the number of factors p in (|4.23p is 
obviously zero. Thus assume that the induction hypothesis holds for a given 1. Then 



[ ^{-^YKbk) W {p + pbpbp)) 2(l+1) ^2(-l) a k(bk) 2a 

(3 a 

= ^2(-lfk(bk) 2l3 (p + pbpbp)j 2 x 

x ^^{-ifkibkf^p + pbpb P )\ 21 Y,(-i) a k(bk) 2a 

13 a 

= ^ — I) 13 k(bk) 213 (p + pbpbp) x (terms with j^p even) 



= (-l) Pl+l32 \(kb) 2f3l p{bk) 2f32 + (kb) 2Pl p{bk) 2f32 bpbp + {kb) 2(3l kbpbk(bk) 2P2 + 
Pi A 

+ (kb) 2 ^ l kbpbk(bk) 2 ^ 2 bpbp x (terms with #p even). 

Now it can be easily checked that when a term with an even number of factors p 
is multiplied from the left by one of the terms in the last sum, then the number of 
factors p either does not change or is raised by two. So in any case one obtains a term 
with an even number of factors p. 

(ii) We again proceed by induction in I. For 1 = 0, there are obviously no contributions 
for p > 0. In the case p = 0, the term of order b 2r in (|4.23p is (— l) r k(bk) 2r . It appears 
exactly 

'r + l — p\ 



1 



times in (|4,23p and has the sign 

(_iy = (-i) 3r +! = 



l 



,4r 2 +3r+ 



l-p 



1 = 



2r + l , 
1 ' 



("I) 2 



This proves the claim in the case 1 = 0. 

Now assume that the induction hypothesis holds for a given I. Consider an operator 
product of the form (|4.24p which involves 2p factors of p. Our task is to determine 
how often and with which sign this operator product appears in the series 



[^i-if k{bk) 2i3 { P + pbpbp)^ (l+1) ^2(-i) a k(bk) 2a 



(4.25) 



Let us first count how often it appears. Since the operator product (|4.24j) is of the 
order b 2r , we only need to calculate up to order 0(6 2r+2 ), 



@25D= [j2(- 1 ) Pk ( bk ? P (P + P b Pbp)) J2(-l) a k(bk) 2a + 0(b 2r+2 ) 
13=0 

r r „, r 

- £■■■ * E- £- + °(' 2 



2(1+1) 

a=0 



(3=0 



(3=0 



a=0 
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r-1 



: ^(-if(kb) 2 ?k + ^2(-if(kb)^kbpb P ] 2 x [ . . . ] 21 ■ ■ ■ +o(b 2r+2 ) 

(3=0 (3=0 (3=0 o=0 
"> ' v ' >< y 



:A =:B 



AA + AB + BA + BB 



L + 0(6 2r+2 ) , (4.26) 



where 



AA=Y, ^(-l) ft+ft (£;&) 2ft p(^) 2/32 , 

(3i =0/3 2 =0 
r r— 1 

AB= ^ ^(-l) /3l+& (A;6) 2/3l p(6A:) 2/32 6p6p, 
/3i =0/3 2 =0 
r— 1 r 

SA= J2(-lf^ 2 (kbf^kbpbkibk) 2 ^ , 

(3 1= 0(3 2 =0 
r—l r—l 

BB=Y Y{-lf 1+P2 {kb) Wl kbpbk{bk) W2 bpbp , 

/3 1= 0/3 2 =0 

r 2 , r 

L= ^{-lfkibkf^ip + pbpbp)^ ^2(-l) a k(bk) 2a . 

(3=0 a=0 

We now treat the three cases p = I + 1 < r, p = 0, and < p < I + 1 separately. 

(a) p = I + 1 < r: 

By the induction hypothesis, the operator product (|4.24|) appears in L exactly 
( r +^) = C^ 1 ) = times. Thus the term (14241) arises in (I4.26P if a summand 
of L which involves 2p — 2 factors of p is multiplied from the left by a term of 
the series E/j=o • • -] 2 ^ n sucri a wa y that two additional factors p are created. 
This can happen in the following ways: 

(1) (^4^4-term ending by k) x (term in L beginning by k) 

(2) (Ai?-term) x (term in L beginning arbitrarily) 

(3) (i?v4-term) x (term in L beginning by k) 

(4) (i?B-term) x (term in L beginning arbitrarily). 

Each term of the form (|4.24p arises as such a product precisely once, because 
in case (1) the first p is at an odd and the second p is at an odd position, 
in case (2) the first p is at an odd and the second p is at an even position, 
in case (3) the first p is at an even and the second p is at an odd position, 
in case (4) the first p is at an even and the second p is at an even position. 
Since by the induction hypothesis, a term involving 2p — 2 factors of p appears 
in L exactly (q) = 1 times, the term (|4.24[) appears in (|4.26p exactly 

I=( r ) = ( r ^ l \ 1) - P ) times. 

(b) p = 0: 

In this case, the operator product (|4.24p is of the form k{bk) 2r . Only an AA- 
term in (|4.26p with 02 = does not create an additional factor p when it is 
multiplied from the right by a term of L involving no factors of p. By induction 



22 



F. FINSTER AND A. GROTZ 



hypothesis, for every r'eN the factor k(bk) 2r ' appears in L exactly ( r t l ) times. 
Thus the factor k(bk) 2r appears in (|4.26[) exactly 



E 



r' + l\ fr + (I + l)-0\ fr + (l + l)-p\ . 
/ ) = [ <l + l)-0 H (/+!)-„ ' 



< p<l + l: 

As in case (a), we want to factor a given operator product (|4.24p at a given 
odd position into a product of a term in E/j=o . • -} 2 and a term in L. Exactly 
as in case (a) one verifies that if this factorization is possible, it is unique. 
Furthermore, this factorization is possible at the first position, at the third 
position, . . . , until the second factor p appears (more precisely, until we are at 
the position of the second factor p, or else the second factor p is at an even 
position and we are at the subsequent odd position). We denote the number 
of such possible factoring positions by A + 1. For 1 < r' < A, factoring at 
the r'-th position yields a term in L of order 5 2 ( r + 1_r ) involving 2p factors 
of p. According to the induction hypothesis, this term in L appears with the 
combinatorial factor (^ r+1 7-p + ^ P ) • Factoring at the last position r' = A + 1 
yields a term in L of order 6 2 ( r_A ) involving 2p — 2 factors of p. According 
to the induction hypothesis, this term in L appears with the combinatorial 
factor ( w^Zi)~ Therefore, the operator product (|4.24|) appears in the 
series (I4.26P exactly 



E 



V + 1 - r') + I - p\ /(r - A) + 1 - (p - 1)\ = fr + (I + 1) - p 
l-p ) + \ l-(p-l) )~\ (l + l)- P 



times. 



We have thus proved the first part of Lemma 14.41 (ii) . 

It remains to prove the induction step for the claim concerning the signs of the terms 
(|4.24p in (|4,26p . By induction hypothesis, an operator product TL of the form (|4,24p 
appears in L with the sign 

(-1) 1+S \ 

where the sum E goes over all positions in TL where Cj = k. This is of course equal to 

/ 2 ^ l+#(factors k at odd positions in H) 

Likewise, each operator product Q appearing in (^4 + B) 2 has the sign 

^ -Q#(f ac tors f. a ^ Qfjj positions in Q) 

When multiplying Q by 7i, a position in the resulting operator product QTL is odd iff 
the corresponding position in the old product Q or TL was odd. We now distinguish 
two cases. 

(a) The term Q ends with a factor k and the term TL begins with a factor k: 

As these two factors k multiply to a factor p by our multiplication rules, the 
resulting term QTL has two factors k less at odd positions than the initial factors 
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Q and H together. Thus the sign of QTt is 

^ _ ^ # (factors k at odd positions in Q) ^_ j^l+#(factors k at odd positions in Ti) 

^ ^^#(factors k at odd positions in 6)+l+#(factors k at odd positions in H)—2 

^ ^^l+#(factors k at odd positions in QTi.) 

(b) In all other cases, one has 

#(k at odd positions in QTt) = at odd positions in Q) 

+ #(& at odd positions in TC), 

and therefore QTi has the sign 

^ -j^#(factors k at odd positions in Q) ^_ ^^l+#(factors k at odd positions in Ti) 

^ 2 ^ l+#(factors k at odd positions in QTi.) 

This concludes the proof. □ 

5. The Unitary Perturbation Flow 

In this section we prove that there exists a operator U, uniquely defined by a per- 
turbation series, which transforms the vacuum operators p m , k m into the interacting 
operators p^ es , k™ s , i.e. 

?T = UpmU- 1 , C S = C/A^- 1 . (5.1) 

The operator U will be unitary with respect to the indefinite inner product on the 
wave functions 

/ '^(x)$(x) d A x , (5.2) 



P 



where ^ = is the usual adjoint spinor. For the proof, we shall consider the 

one-parameter family of rescaled operators 

r5,:=(pD T ^d fc-=(C) T (5-3) 

corresponding to the family of external fields (t25) t >o, i.e. 

(i$ x + rS(x) - m) ^(x, y) = = + rS(x) - m) y) . (5.4) 

The idea is to consider the differential equations for 

the so-called perturbation flow equations, and to write them purely in terms of in- 
teracting operators. This will reveal that the perturbation flow equations have the 
commutator structure 



d 



A{T).f>l, and ^.Kn= Mt).K, ( r ».5) 



with an anti-symmetric operator A, from which U can be obtained by integration. 
In preparation, we define the operator 

noting that this implies the relations 

s±=s m ± mp m . (5.7) 
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Moreover, we can generalize the formulas (|2.15p and (|2.16p . 



Lemma 5.1. The operators j5^ s , k™ s and s m satisfy the following identities: 



PP 

^res ~ ~ -res _ fi™s 

Pm a m' — »m' Pm ~ i Pm 

m — ml 

PP _ 

J, res ~ _ ~ , jures _ Lres 

^m °m' — °m' ^m ~ _ , ^m 



Proof. In view of Theorem 14.14 it suffices to show that 



> ~ PP > < PP < 

Pmb m s~m' = 7 Pm,b m and s m i b m p m = b p r 

m — m m — ml 



Substituting the relations 



&m — i i^Pm 



into the power series (I4.12P of s m , we note that the terms with an odd number 
factors p m cancel in (|5.6p . We thus obtain 



5 m =^(-l) n (™)* Q C m (Q,l)'BC m (Q,2)...'BC m (Q,n + l) 

n>0 QeV(n+l) 
#Q even 

= E(- 1 )™ E (™)* Q s m HC m (Q,2)...T>C m (Q,n + l) 

n>0 Q£V(n) 
#Q even 

+ E(" 1 ) n E (^)* Q Pm 'BC m (Q,2)...T>C m (Q,n + l), 

n>0 Q&V(n) 
#Q odd 



where V(n + 1) is defined as in Theorem 14.11 and 

C m (Q,j) : 



p m , j e Q 

s m i j ^ Q 



Moreover, we introduce the abbreviations 



E : = B- 1 )" E and E : = E(- 1 ) n E 

n,even n>0 Q£V(n) n,odd n>0 QeV(n) 

#Q even #Q odd 



THE RESCALED CAUSAL PERTURBATION EXPANSION 25 

Applying the calculation rules from Lemma 12. 11 we obtain 



k>{) 



■ Pm Sm'3C m /(Q,2)...3C m /(Q,n + l) 

n,even 

fc>l n,even 

+ PmJ] p m '3C m /(Q,2)...3C m /(Q,n + l) 

n,odd 

-Pm^(-3s m ) fc - 1 3s m p m >'BC m ,(Q,2)...'BC m ,(Q,n + l) 

k>l n,odd 

' >T> 7 p m SC m ,(Q ) 2)...3C m /(Q,n + l) 



m — m 

n, even 



PP 

yPm^-^m)^™ X] 3C7 m ,(Q,2)...SC m /(Q,n + l) 



m — m 

k>0 



PP 

+ 7 Pm V(-£s m ) fc Bs m , V BC m /(Q,2)...BC m /(Q,n+l) 

k>0 n,even 



- vr 2 5(m - ml) p m ^(-3s m ) fc S ^ p m 2C m (Q, 2) . . . £C m (Q, n + 1) 

fc>0 n,even 

+ <y(m - ml) Pm%C m (Q, 2)... T>C m (Q, n + 1) 

i,odd 

/3m ^(-£ Sm ) fe £ £ p m ^C w /(Q,2)...3C m /(Q,n + l) . 



n,odd 

PP 



ml — m 

k>0 n,odd 

Now a careful inspection shows that in the last equation, the fourth and the fifth line 
cancel each other, whereas the k = 0-terms in the third and the sixth line cancel the 
n > 1-terms in the first line. Furthermore, the k > 1-terms in the third and sixth line 
cancel the n > 1-terms in the second line. We thus conclude 

pp pp 

Pm&m 8m' = l(Pm + Pm V"(-Ss m ) fe ) = lPmb m ■ 

m — ml m — m 

k>l 

The equation s m > b^Pm = ^ m can be proven analogously. □ 

Next we need to derive a completeness result for the spectral projectors of the inter- 
acting Dirac operator. The completeness relation could be stated as 



p™ s dm = 1 . 

However, since we introduced the spectral projectors only for a real and positive mass 
parameter m, it is more convenient to work instead with contour integrals of the 
Green's function in the complex plane. To this end, we introduce for any e > the 
hyperbolas 

C £ = {n G C with |(Re/i)(Im/i)| = e 2 } . 
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We also consider C e as a contour along which we integrate in anti-clockwise orientation. 
In view of the identity (p 1 ) 2 = p 2 £ R, the eigenvalues of the matrix p 1 lie in R U iR. 
Hence for any /j, £ C e , the free Green's function can be defined in momentum space 
by 

^(p) = (^-/x)- 1 . (5.8) 
Moreover, the Green's function with interaction can be defined by the perturbation 
series 

oo 
n=0 

because writing the Feyman diagrams in momentum space, one gets products of the 
form 

S M (p n ) ^(Pn ~ Pn-l) • ' • S M (jPl) S(pi - Po) ^(Po) (5-10) 

where each factor is well-defined according to (|5.8p . Next, we define the principal value 

by 

X- M-s ' " ' ^ ~ Ic- 2 (/U-a; + 4ee^ + /z - x - 4e ) " ' Sm ^ ' 

where the phase </3 = y?(x) is to be chosen such that the points x =t Aee lLp are both 
outside the set 

J £ :={/i£C with |(Re/i)(Im/z)| < e 2 } 

(in the limit e \ 0, the results of all contour integrals will become independent of the 
choice of (p). 

Lemma 5.2. For any [i £ C £ and m > 0, the following identities hold to every order 
in perturbation theory, with convergence in the distributional sense: 

— lim / s M dfx = 1 (5.12) 



2-ni e \o 

lim f 

2iri e\o Jq b fx — m 



lim / — s M dfj, = s m (5.13) 



m — fi m — ix 

Here we used the definitions (|5.8|) . (|5.9|) and (|5.11|) . whereas s m is defined by (|5.6|) . 

Proof. In order to derive (|5.12p . we consider the perturbation series in momentum 
space. Using again that the spectrum of the matrix p 1 lies in R U iK, we can inte- 
grate (|5.8|) with residues to obtain 

— — = -2-iri 1 . 
c e V ~ V 

Similarly, when integrating (|5.10p . the integrand is a product of poles (p — x)~ p of 
total order larger than one, and the poles are all enclosed by the contour C e . Hence 
computing the integral with residues, we obtain zero. This proves (|5.12p . 

For any \i £ C £ and v £ C \ J £ , a straightforward calculation using (|5.8p and (|5.9p 
shows that to every oder in perturbation theory, the "resolvent identity" 

1 . . 
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holds. Integrating \x on both sides over C £ and using that 

s^s u dn = (J dfj\ s u -2m s v 

l r l 

a\x = s u \ OLfi = , 



we find that 



c e H ~ v Jc e p-v 

L dyu = — 2m' s„ . 



Now (|5.13|) follows immediately from the definition of the principal value (|5.1ip . 
Similar as in the proof of Lemma l5.ll the identities (15.141) follow from the relations 

•5 ab rn Pm — b m Pm and s nb m k m — 6 fe m . 

m — fjj m — fj, 

These relations follow immediately from the definition of 6^ , (|2,12p , by applying the re- 
solvent identity (|5.15p (with v replaced by m) as well as the relations (I2.15P and (|2.16p , 
which after replacing m' by the complex parameter fj, £ C E are clearly valid even with- 
out the principal part. □ 

We are now ready to state the equations for the perturbation flow. Using a notation 
similar to (|5.3p and (|5.4p for the Greens functions corresponding to the family of 
external fields (rS) r >o, we write 



A(t).i>1 and ^Kn = A ^)-K, 



Theorem 5.3. The perturbation flow equations have the commutator structure 

dT Pm 

where A(t) is a family of operators, which are anti- symmetric with respect to the inner 
product (|5.2p . 

Proof. Since the Greens functions (^j) T are defined via the formal Neumann series 
(I4.12|) . their r-derivatives are given by 

Using the definition of p T m , (|4.11[) . and applying the relations ()5.7fl . we obtain 

JjlPm = -S T m Wm ~ Pm^Ki ■ ( 5 ' 16 ) 



The identity 



= -^(i0 + r'B-rn)kl l = 'Bkl l + (i$ + rS - rnj ^-k T m 



implies that the r-derivative of k T m is a solution of the inhomogeneous equation 

%$ + r% - ^—k T m = -^k T m . 

Thus, the r-derivative of is of the form 

d_ 7 
~drr 



k T = — 5 r e nk T -a- r H~ r 
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where "H^ is a solution of the homogeneous equation 

'i0 + T , B-m)y(Z n = Q. (5.17) 



Since is a symmetric operator, the same is true for its r-derivative. We thus 
conclude that 

Tf" = — h T 5 T -i- < W" r 

where r K T m is another solution of the homogeneous equation (|5,17p , which is moreover 
symmetric with respect to the inner product (15, 2j) . The most general ansatz for "K m is 

r H" r = fT A T fT + k T C T k T + k T T> T r> T 4- f> T (T) T Yk T 

J L m fm^mfm ' n 'myrn n "m ' "'m rafra T ™r mi ""m 

with a certain operator and certain symmetric operators _4.J„ and C^. This leads 
to the identity 

— (h T h T \ — —h T r R« r £ T — £- T 5 T r R£- T 



-\-5(m m ) ( s m r Rp m p m 'Bs rn + p m A m k m + k m C m p m + k m T^ m k m 

4- r> T ('T) T ^rj 1 " -I- t r i T r> T -I- rf f r fr T -I- rT T) t r> T 4- £ T fT) 1 " 
^/'mV^m/ ^ ""m^mrm ^ fmym n "m ' I'm'-' mfm ' n "m\ 1 -' m) n . 

In view of Lemma [5.1l the two terms in the first line cancel each other, and the identity 

d ~ ~ d 

-j-{k T m k T m ,) = 5(m - m')—p T m = 6(m - m'){-s T m 'Bp T m - p T m ^~s m ) 
yields the relations 

A m = —C m and v^m) = • 

Thus 

~7^m = ~ ^rn^^m ~ ^rn^^m Prn^mPm ~ ^m-^m^m + ^mPmPm ~ PmPm^rn (5.18) 

with 

(A m )*=A m and ra* = -2^. (5.19) 
We now define the operators 

i r r pp 

M T ) '■= — o lim / du I dv slUsT. 

47T 2 e\0j Cs P J Ce V-H " 

J / (fi T T) T n T — h T T) T h T 4- h T A T f) T — f) T A T h T I dm 

2 J yrm'-^rnl'rn n rn 1 ^rn n "in ~ ^rn^rnfm Pm- n m n 'm J Ujl "' ■ 

From (|5.19p it is obvious that A{r) is an anti-symmetric operator. Furthermore, ap- 
plying Lemma 15.21 we find that 



PP f f PP 1 



dfj, / du s Tr Bsl p T m = dfi du s Tr B p T m 

c e Jc E v — J Ce J Ce v- fj, * m-v 

f PP 

= -2m / dfi s Tr Bp T m = -^ 2 s m <Zp m ■ 

Using this identity together with the analogous identity with p T m replaced by k m , one 
readily verifies that (|5. 16[) and (|5.18h can be written as 

= M T )Pm + Pm A ( T )* and ^k T m = A(r)k T m + k T m A i T T , 
concluding the proof. □ 
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As mentioned at the beginning of this section, the result of Theorem 15 . 3 1 can be used 
to prove the existence of a one-parameter family of unitary transformations describing 
the interaction. 

Theorem 5.4. There exists a one-parameter family of operators U(t), which are 
unitary with respect to the inner product (|5.2p , such that for any mass m > and for 
allT>0 

P T m = U(T)p m U(ry and k T m = U(r)k m U(T)* . (5.20) 
Proof. If there exists such a family U(t), then taking the T-derivative of (|5.20p yields 

P(r)p m (r)(/(r)- 1 - U{r) Pm U{T)- 2 U{r) 



On the other hand, from Theorem 15.31 we know that the perturbation flow equations 
do have the commutator structure f)5.5f) . Thus integrating these equations and taking 
ordered exponentials yields a family of unitary operators 

U(t) = Texp ( - J A(t) dr) (5.21) 

fulfilling the identities flOOj) . □ 

We remark that the operators A^ and T> T m are uniquely defined by perturbation 
series, for a perturbation being of the order ©CB 1 ). The same is then true for the 
operators A(t). This implies that the operator U(t) in (|5.2ip has a well-defined per- 
turbation expansion. Going through the combinatorial details, it would be straight- 
forward to derive explicit formulas for this perturbation expansion to any order. Here 
we are content with the existence statement which immediately implies the following 
interesting corollary. 

Corollary 5.5. There exists a unitary operator U, defined by a unique perturbation 
series, which for any mass m > transforms the free fermionic projector (|1.3|) into 
the interacting fermionic projector ()4.ip : 

P m = UP^ a U* . (5.22) 

6. The Light-Cone Expansion of the Fermionic Projector 

The light-cone expansion is a useful technique for analyzing the fermionic projector 
near the light cone. An introduction to this technique can be found in [9], for details 
see [U [8] . In this section we will show that our rescaling procedure has no effect on 
the singularities of the fermionic projector on the light cone, but it does change the 
regular so-called high-energy contribution. 

We briefly recall the definition of the light-cone expansion. 

Definition 6.1. A distribution D(x, y) on K 4 x R 4 is of the order 0((y — x) 2p ) on the 
light cone, p£Z, if the product 

(y-x)- 2 PD(x,y) 



d 



dr m dr 



and 



dr m 
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is a regular distribution, i.e. a locally integrable function. It has the light cone expansion 

oo 

D(x,y) = Y,D [j] (x,y) (6.1) 

3=9 

with g £ 7L, if the distributions D^\x,y) are of the order 0{(y — x) 2j ) on the light cone 
and if D is approximated by the partial sums in the sense that 

p 

D(x,y)-Y,D [ J ] (x,y) (6.2) 

3=9 

is of order 0((y — x) 2p+2 ) on the light cone for all p> g. 

We note that the lowest summand D^(x,y) gives the leading order on the light cone. 
If D is singular on the light cone, g is negative. 

By explicitly expanding the series (|2.8p term by term and carrying out the sum over 
the orders in perturbation theory, in [8] the light-cone expansion of the causal Green's 
functions s v and s A is performed. Using (|2.2ip . one readily obtains the light-cone 
expansion of k. The so-called "residual argument" [8_, Chapter 3] allows us to carry 
over the light-cone expansion to the operator p ms , (|4,2p . In this way, in [8] the light- 
cone expansion is derived for the so-called residual fermionic projector P res defined 
by 

P^:=l(p^-~k). (6.3) 

The residual fermionic projector is a solution of the Dirac equation (jl.ip . but its pertur- 
bation expansion has a different combinatorics than that of the fermionic projector P 
as given by (|4.ip . The difference of these two operators is defined to be the high-energy 
contribution P he (x,y), 

phc. = p_pres (g 4 ) 

The following proposition extends the result of [8J Theorem 3.4] to the rescaled setting. 

Proposition 6.1. Under the assumptions of Lemma \2.6[ the integral kernel P he (x,y) 
of the high-energy contribution is, to every order in perturbation theory, a smooth 
function in x and y. 

Proof. Using the formulas (14. ip and (16. 3p for P and P res ; one sees that 

P he = \(k-k* es ) . 

As will be shown below, the following replacement rule holds: 

If in the perturbation series (|4.3p for k Tes one replaces all factors p by 
factors k, one obtains precisely the perturbation series (|2.22p for k. 

Furthermore, we know from Lemma 14.31 that the total number of factors p in every 
summand of the perturbation expansion for k Tes is even. Hence to every order in 
perturbation theory, k can be obtained from k res by iteratively replacing pairs of 
factors p by factors k. Thus, to every order in perturbation theory, the operator P hc 
can be written as a finite sum of operator products of the form 

C^B • • • C CT+1 :B (k-BC^x ■ ■ ■ Cr+xKk - pKC^ ■ ■ ■ C T+ iSp)sC T _i ■■■'BCo. (6.5) 
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It is shown in the proof of [H Theorem 3.4] that terms of the form (|6.5p are smooth 
and bounded in position space, provided that the assumptions of Lemma 12.61 hold. 

It remains to prove the above replacement rule. For convenience, the replacement 
of all factors p by factors k is denoted by the symbol ~~>. Furthermore, we introduce 
the abbreviation 

S(Z) : = [ Y,(- i7T ) 2(3k (bk) W {p + TT 2 P bpbp)\ 21 J2(-^) 2a Hbk) 2a . (6.6) 
In the proof of Lemma 14.31 we saw after (|4.20p that for any r 6 N, 

r 

~k™ = J2 fl,r b K §(0 b> + 0(6 2r+2 ) . (6.7) 

1=0 

Replacing in (|6.6p all factors p by k, we obtain 

S(Z) ~» [ ^(-ivr) 2/3 (A;6) 2/3 (l + vr 2 ^)] ^(-iit) 2a k{bk) 2a + 0(& 2r+2 ) 

/3=0 a=0 
r 

= ^{-iir) 2a k{bk) 2a + 0(6 2r+2 ) , 

o=0 

where in the last step we used that the sum in the square bracket is telescopic. Sub- 
stituting this formula into (|6.7p . we obtain 

r r r 

~k Tes - £ fi,r £Htt) 2q &< fc(6fc) 2a b> + 0(6 2r+2 ) = £ / l)P fc + 0(6 2r+2 ) , 

Z=0 Q=0 «=0 

where we used (|2.22|) . The result now follows from the identity 

r 

£Ar = l, 

1=0 

which is verified by an elementary combinatorial argument. □ 

According to (|6.4|) and the above Proposition, the singular behavior of the fermionic 
projector P is completely described by the residual fermionic projector p rcs . In [5] 
the residual fermionic projector is decomposed into two parts, 

Thus, the light-cone expansion of the fermionic projector becomes 

P(x, y) = P causal (x, y) + P lc (x, y) + P he (x, y) . 

Here the contribution p causal encodes the singular behavior of P on the light cone, 
whereas P he (x,y) and P le (x,y) are, to every order in perturbation theory, smooth 
functions in x and y. 

We finally explain how the "causality" of the causal perturbation expansion is to 
be understood. Causality entered our construction when defining k in terms of the 
causal Green's functions (12.21[) . This is reflected in the kernel P causal (3;, y), which is 
local and causal in the sense that it is given explicitly in terms of integrals of H> and 
its partial derivatives along the line segment xy. However, taking the absolute value 
of an operator (|2.25p as well as the rescaling procedure (|3.5p do not preserve causal- 
ity, but these operations are necessary in order to distinguish a properly normalized 
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subspace of the Dirac solution space. The non-causality of our construction becomes 
apparent in the high- and low-energy contributions, which involve integrals of B over 
the whole space-time. Thus our construction reveals that, although the Dirac equation 
is causal, the interacting Dirac sea is a global object in space-time, violating locality 
and causality. For a discussion of this remarkable fact we refer the reader to the book 



Appendix A. The Leading Orders of the Perturbation Expansions 
The operator k in (|2.22p is given by 

k =k - sBk - kBs + kBsBs + sBkBs + sBsBk - i^kBkBk 
- kBsBsBs - sBkBsBs - sBsBkBs - sBsBsBk 

+ 7T 2 (sBkBkBk + kBsBkBk + kBkBsBk + kBkBkBs) + OCB 4 ). 



The operator X in (|2.44j) has the expansion 

X =vr 2 [p'Bp'Bp - pT>k-Bk - k'Bk'Bp + k'Bp'Bk 

+ p'Bs'Bk'Bk - k'Bp'Bs'Bk - p'Bp'Bs'Bp + p'BkSsSk 

+ k'Bk'Bs'Bp - k'Bs'Bp'Bk + k'Bs'Bk'Bp - p'Bs'Bp'Bp) + 0('B 4 ), 



p=p- sT>p - p"Bs + pBsfts + s'Bp'Bs + s'Bs'Bp 

+ y (k'Bp'Bk - kBkBp - pBkBk - pBpBpj 

- pBsBsBs - sBpBsBs - sBsBpBs - sBsBsBp 
vr 2 / 

+ — [sBpBpBp + pBsBpBp + pBpBsBp + pBpBpBs 

+ kBsBkBp + kBkBsBp + pBkBsBk + pBsBkBk - kBsBpBk 

- kBpBkBs - kBpBsBk + pBkBkBs + sBpBkBk + sBkBkBp 

+ kBkBpBs - sBkBpBk) + 0(B 4 ) , 



Y =p - vr 2 [pBpBp - pBkBk - kBkBp + kBpBk 

+ pBsBkBk - kBpBsBk - pBpBsBp + pBkBsBk 

+ kBkBsBp - kBsBpBk + kBsBkBp - pBsBpBp) + 0(B 4 ). 



The operator p Tes in (|4.10p has the expansion 

p vcs = p - sBp - pBs + pBsBs + sBpBs + sBsBp - i^pBpBp 
- pBsBsBs - sBpBsBs - sBsBpBs - sBsBsBp 

+ vr 2 (sBpBpBp + pBsBpBp + pBpBsBp + pBpBpBs) + 0{B 4 ), 



which yields for the operator p in (|2.39|) 



and for the rescaling operator Y in f|3.3[) 
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whereas the operator k TCS in (|4,16p is given by 

k ies =k - sBk - kBs + kBsBs + sBkBs + sBsBk 

IT 2 ( 

+ ^P^P + P'Bk'Bp - pBpBk - kBkBk 

- kBsBsBs - sBkBsBs - sBsBkBs - sBsBsBk 

7T 2 / 

+ — [sBkBkBk + kBsBkBk + kBkBsBk + kBkBkBs 

+ kBsBpBp - sBpBkBp + sBpBpBk + pBsBpBk - pBkBpBs 
+ kBpBsBp + kBpBpBs - pBkBsBp + pBpBkBs + sBkBpBp 

- pBsBkBp + pBpBsBk^j + 0(B A ). 
Thus for the operator t in (|3.ip we obtain 

1 



< 2 



p-k - sBp - pBs + sBk + kBs 
+ pSsSs + sBpSs + sBsBp - kBsBs - sBkBs - sBsBk 



+ -K 2 kBkBk + y (kBpBk - kBkBp - pBkBk - pBpBp 

- pBsBsBs - sBpBsBs - sBsBpBs - sBsBsBp 
+ kBsBsBs + sBkBsBs + sBsBfcSs + sBsSsSA; 

2 

+ — (sBpBpBp + pBsBpBp + pBpBsBp + pBpBpBsJ 

- ir 2 (sBkBkBk + kBsBkBk + kBkBsBk + kBkBkBs^j 
vr 2 / 

+ y (fc'Bs'B/cBp + fcB/fc'Bs'Bp + pBkBsBk + pBsBkBk 

- kBsBpBk - kBpBkBs - kBpBsBk + pBkBkBs 

+ sBpBkBk + sBkBkBp + kBkBpBs - sBkBpBk) 
and the rescaled fermionic projector (|4.ip is given by 
p-k - sBp - pBs + sBA; + kBs 



+ 0(B 4 ), 



P = i 
2 



+ pBsBs + sBpBs + sBsBp - kBsBs - sBkBs - sBsBk 

IT 2 / 

- 7r 2 pBpBp + — [kBkBk + pBpBk - pBkBp + kBpBp 

- pBsBsBs - sBpBsBs - sBsBpBs - sBsBsBp 
+ kBsBsBs + sBkBsBs + sBsBkBs + sBsBsBk 

+ vr 2 (sBpBpBp + pSsSpSp + pSpSsSp + pBp'Bp'Bs 

- y (sBkBkBk + kBsBkBk + kBkBsBk + fcBfcB&lBs) 
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7T 2 



+ — [p'Bk'Bp'Bs - s'Bk'Bp'Bp + p'Bk'Bs'Bp - p'Bp'Bk'Bs 
- p'Bs'Bp'Bk + p'Bs'Bk'Bp - k'Bp'Bs'Bp - k'Bs'Bp'Bp 



p'Bp'Bs'Bk + s'Bp'Bk'Bp - k'Bp'Bp'Bs - s'Bp'Bp'Bk 



+ OCB 4 
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